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Through a determination of the 5-matrix element for the scattering of a localized 
colour dipole on a proton, we show that saturation effects can already be sizable 
at HERA. The saturation scale is found to be around 1-1.5 GeV^ for the highest 
available energies and for central collisions, and is a decreasing function of the 
impact parameter. 



Presently, there are basically two most discussed theoretical descriptions 
of the HERA data for the proton structure functions, equally successful at 
small- a; but contradictory in their physical implications. On one hand, models 
based on the DGLAP equation fit the data very well. On the other hand, the 
recent saturation modelo reproduces all small- a; data with a small number 
of parameters. The latter assumes a unitarized cross section but no scale 
evolution is taken into account. Its success could mean that saturation effects 
are already important at HERA, at a scale around 1 GeV^. However, the 
very question whether we have reached or not the regime of high densities at 
HERA in the perturbative domain has so far not been addressed in a direct 
way. We propose here a simple method to determine reliably how dense the 
proton looks. 

1 A brief introduction to saturation 

The HERA physics in the kinematical range considered here is the physics of 
photon-proton interactions. We will denote W the 7*p center-of-mass energy, 

the virtuality of the photon and x ~ /W^ the Bjorken variable. In a 
Bjorken frame in which the proton has a large light-cone momentum 0, 0) 
and the photon the momentum (0, , the photon-proton interaction can 
be interpreted as the photon probing the partonic content of the proton. This 
probe is sensitive to collinear fluctuations down to transverse distances of order 
l/q-*- oc l/Q, and resolves fluctuations occuring over times as short as l/q~ oc 
p'^ /W'^ . Thus by increasing Q^, one reveals more collinear partons. By raising 
the center of mass energy W , one becomes sensitive to softer fluctuations. 

If one keeps increasing W at fixed the density of partons of size 1/Q 
will eventually be so large that their wave functions will overlap, and recombi- 
nation effects will become sizable. The borderline between these two regimes 
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is usually parametrized by the function = QKW"^). Qs is called the "sat- 
uration scale" . The question is whether we have already crossed this line at 
HERA. 

Generally speaking, wave diffraction allows to obtain a picture of a micro- 
scopic object. A mere Fourier transform then relates the diffractive pattern to 
the density profile of the object. This can be done at HERA using diffractive 
processes, as will be discussed in the following. 

2 Transparency of the proton to dipoles from difFraction 

At high energy, one can relate diffractive vector meson electroprodbiction to 
dipole-proton elastic scattering. This is the celebrated dipole modelH, that we 
briefly recall here. 

In the proton rest frame, the scattering proceeds in three steps. The 
virtual photon first breaks up into a qq pair of size r, which then scatters off 
the proton, before being recombined into a meson. High W guarantees that 
these processes are well separated in time. At HERA, the qq pair travels over 
distances of a few 10^ fermis before entering the interaction region, of size 1 
fm. The recombination takes a comparable "time" . This remark leads to the 
following factorization formula for the amplitude A of this process: 

A{Q',W',A) = I d^rA{W\r,A) J Vq(^, r)Vy(2, r) , (1) 

where A is the dipole-proton elastic amplitude, and are the photon and 
vector meson wave functions respectively. The latter has to be modelled: we 
considered several models available in the litterature. A is the 2-momentum 
transfer: it is conjugated to the impact parameter, and this will be the variable 
with respect to which we will take the Fourier transform. We define t = |A|^. 

Thus the proton sees a beem of asymptotic qq dipoles, distributed ac- 
cording to the overlap of the photon and meson wave functions (second in- 
tegral in formula ([l])). On very general grounds, one can argue that this 
function is sharply peaked around a value rg of the dipole size scaling like 
I/a/Q^ + My {My is the mass of the meson, see El). Hence the amplitude 
can be approximated by ^(Q^, A) ~ A{W^,rQ,A) x N{Q) where N{Q) 
is the normalization of the dipole distribution. Note that N{Q) is a priori 
model-dependent, but in practice this dependence is very weak. 

Once we have extracted the dipole-proton amplitude from the HERA data 
in this way, we perform a Fourier analysis, by defining the function 

S{W',rQ,b) = l-l J ^e-^^^AiW^rQ,A). (2) 
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One easily checks that S{W'^, r, b) is the S'-matrix element for the scattering 
of a localized dipole of given size r at given impact parameter b. This holds in 
particular because at high energy such a dipole state can only be absorbed: the 
^-matrix is diagonal and its eigenvalues are real between and 1. They can 
be interpreted as kind of "transparency coefficients" . Furthermore, note that 
l — S'^{W^,r,b) is the probabihty for such a dipole state to undergo inelastic 
scattering. 

3 A picture of the proton at HERA and its interpretation 

We have applied the procedure described above to the HERA data. We took 
the available data for diffractive production of longitudinal p mesons from the 
HI collaborations, x was taken small, around 5 • 10~^. Three different values 
of were considered, = 7, 3.5, and 0.45 GeV^, corresponding to dipole 
sizes of about 0.16, 0.21 and 0.35 fm respectively. The data was available only 
for t < 0.6 GeV^, so we had to extrapolate them to larger values of A before 
performing the Fourier transform. We made three such assumptions. The 
results are represented on fig. |l|. We get a reliable extraction of 5(6) down to 
b = 0.3 fm. 

We see that for say < 2 GeV^ and b < 0.4 fm, the probabihty of 
inelastic interaction is more than 50%. Hence the proton looks already quite 
black in its central region at present HERA energies. To see whether this 
"blackening" occurs in the perturbative regime, we estimate the saturation 
scale. For this purpose, it is necessary to choose a parametrisation for S. Let 
us consider the following Glauber-like form: 

SiW^r,b)^e-Q''^'^"-'>>"^\ (3) 

where the fact that Qs controls the onset of the saturation regime is quite 
explicit. From this formula we now extract the saturation scale and obtain 
Q2 ~ 1 - 1.5 GeV^ for b = 0.3 fm ancLQ^ ~ 0.2 GeV^ for b = 1.0 fm. 

We refer the reader to the paper □ for all the details of this analysis and 
for more results. 
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Figure 1. S-matrix element as a function of fe, for x = 5 ■ 10~*. Three values of 
are considered, and three different extrapolations of the data at large momentum transfer 
{exp(— At), The sensitivity to this extrapolation is high in the grey zone. 
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